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Abstract. The dynamics of gaseous stars can be described by the Euler- 
Poisson system. Inspired by Rein's stability result for 7 > ^, we prove the 
nonlinear instability of steady states for the adiabatic exponent 7 = g in 
spherically symmetric and isentropic motion. 



1. Introduction 

The motion of gaseous stars can be described by the Euler-Poisson equations: 

(1.1) | + V'H = o 

(1.2) + y • (pu ® u) + Vp = -pV$ 

at 

(1.3) A$ = 47175 

where t > 0, x E R 3 , p is the density, u E R 3 the velocity, p the pressure of the 
gas, and <& the potential function of the self-gravitational force. We consider the 
isentropic motion i.e. p — Ap 1 ^ < 7 < 2), where A is an entropy constant and 
7 is an adiabatic exponent. In our case A will be normalized as =2- and 7 will be 

chosen as f . 

5 

For the spherically symmetric motion the above equations where p > can be 
written as following: 

ii dp dp du 2 

(1.4) -£ +u ^ +p + _ pu = , 

at or or r 

. du du I dp 47T f r o , 

1.5 —+u — + - 7 f + — p S 2 ds = Q. 

01 or p or r z J 

First consider the stationary solutions (po(r), 0) of (1.4) and (1.5). They satisfy 

(1 „ i + t?£s*-o. 

This ordinary differential equation has been well studied. One interesting ques- 
tion relevant to our context can be given like this: Given the mass M > 0, how many 
solutions are there for (1.6) with J R3 pdx = M(p) = Ml For our purpose we sum- 
marize the answer according to the range of 7. See [3] or [7] for details. If | < 7 < 2 
and any M > 0, there exists at least one compactly supported stationary solution 
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p such that M{p) = M. For | < 7 < 2, every stationary solution is compactly 
supported and unique. If 7 = | and any M > 0, there is a unique ground-state 
solution (not compactly supported) p. The solution can be expressed in terms of 

Lane-Emden function and moreover it can be written explicitly as On (r) = — 5- 

(l+r 2 )3 

up to scaling in r and constant multiplication. On the other hand, if 1 < 7 < |, 
there is no stationary solution with finite total mass. 

The stability question has been a great interest and it has been conjectured by 
astrophysicists that stationary solutions for 7 < I are unstable; indeed one can 



3 

easily check that when 7 <E (1, |) steady states are not minimizers of the energy 
functional 

1 „nj„ 1 f f p( x )p(y) 
y\ 

by constructing a scaling invariant family of steady states, and this indicates the 
possibility of certain kind of instability. So far only partial results are known in 
this direction. The linear stability of the above stationary solutions was studied 
in in Lagrangian formulation. It was shown that any stationary solution is 
stable when 7 £ (|, 2) and unstable when 7 £ (1, |). In accordance with the linear 
stability, a nonlinear stability for 7 > | was established recently in |S] by Rein using 
the variational approach based on the fact the steady states are minimizers of the 
energy functional E defined in (1.7). See jHj for a great overview of mathematical 
results on the nonlinear stability problems. For 7 = |, the energy of a steady state 
is zero and any small perturbation can make the energy positive and cause part of 
the system go off to infinity which implies an instability of such a state. This kind 
of instability was investigated in 0]. However, the method they employed is not 
applicable when 7 < |. The stability question for | < 7 < | under the physical 
consideration waits for a satisfactory answer. 

Our main result in this paper concerns a fully nonlinear, dynamical instability 
of the steady profile for 7 = | : 

(1.8) Po (r) = 1 and u (r) = 

(1 + r 2 )^ 

We shall show that this steady profile is unstable under an appropriate energy- 
like measurement Ei which will be precisely defined later in Section 5. Here I < 
represents the strength of weights. As / gets smaller, E\ is equipped with the stronger 
weight. Indeed, £0 corresponds to the positive part of the real energy E. First we 
rewrite the Euler-Posisson system (1.4) and (1.5) for 7 = | by letting p = po + a : 

(1.9) a t + -^(r 2 p u) r + -^{r 2 cru) r = 

(1.10) u t + uu r + ^r(po + o~)~Hpo + er) r + [ {Po + cr)s 2 ds = 

lt> T Jo 

If we define and T s by 

= / (Po + a ) u2 + tt(Po + aY^a 2 dx and 

JR3 15 

1 

T s = —= In - where t/JIq is the sharp linear growth rate, 



(1.7) E = E(p, u) = J (ipM 2 + -^P)dx - \ 



-dxdy 
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the main theorem can be stated as follows : 



Theorem 1.1. Let I < —3 be fixed. There exist 9 > and C > such that for any 
small S > 0, there exists a family of solutions f S (t) = of (1.9) and (1.10) 

such that 

(1.11) V^(0) < CS, but sup v/ij?(f) > 9. 

0<t<T 6 V 

We remark that the escape time T s is determined by the exponential growth 
rate of the linearized system (2.1) and (2.2) and the instability occurs before the 
possible blowup of the smooth solutions. The local existence of regular solutions 
for 7=| including (1.8) was shown in jS] and independently. Gamblin uses the 
paradifferential calculus and Bezard uses the Gagliardo-Nirenberg inequality and 
Littlewood-Paley theory. In this paper we take their existence results for granted 
without proving again. 

Rein's work in [5] implies that the steady states for 7 > § are stable under the 
energy functional E as minimizers. By contrast with those cases, our steady profile 
for 7=| docs not minimize E and that motivated us to investigate an instability 
under the energy- like measurement £q. The proof of Theorem 1.1 is based on the 
bootstrap argument from the linear instability to the nonlinear dynamical model. 
In general, passing a linearized instability to a nonlinear instability needs much 
effort in the PDE context: the spectrum of the linear part is fairly complicated 
and the unboundedness of the nonlinear part usually yields a loss in derivatives. In 
order to get around these difficulties, the careful analysis on the linearized system 
which controls the sharp exponential growth rate of its solutions is necessary. In 
addition, the energy estimates on the whole system and the interplay with the linear 
analysis can close the argument. See [H] for its original method. 

The main difficulty in this paper is to derive Proposition 6.1, a key estimate for 
the bootstrap argument. There are two important ingredients: the first idea is to 
find the form of 

(1.12) < n£ + (£) 2 + (quadratic but lower derivative terms) 

where r\ is smaller than the sharp growth rate ^JJio, but it turns out that (1.12) is 
not enough to close the energy estimate; so we introduce the weighted energy £\ by 
the utilization of a family of symmetrizers as weights. New undesirable quadratic 
terms come out during the weighted energy estimates, but they turn out to have 
weaker weights. As for weighted cubic terms, we introduce the weighted Gagliardo- 
Nirenberg inequality (6.19). Here is the modified version of (1.12): 

^ 5~ + (quadratic but better weighted terms) 

+ (£i) 2 + (quadratic but lower derivative terms) 

The potential part V$ not only has the smoothing effect to the whole system 
but also behaves nicely with respect to weights. Eventually, in cooperating with 
weights, the chain type estimates by using (1.13) complete the bootstrap argument 
as well as the proof of Proposition 6.1. 
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Another technical difficulty lies in that both linearized and full system do not 
comply submissively with spatial derivatives. Our argument in both linear and non- 
linear parts heavily depends on the estimates on pure temporal derivative terms. 
As for spatial and mixed derivative terms we use the equations directly to control 
them in terms of temporal derivative terms. Furthermore, when one uses the polar 
coordinates, we seem to end up with simpler one dimensional flow but the singu- 
larity at the origin comes into play. For clear and precise understanding of the 
largest growing mode, we take the polar coordinates in the linear analysis while 
the rectangular ones are used for the nonlinear analysis to take care of spatial and 
mixed derivative terms. 

Other cases of 7, | < 7 < | attain more physical, interesting feature involving 
the vacuum boundary. For those cases, steady states satisfying (1.6) are compactly 
supported and even the local existence of the Euler-Poisson system including those 
stationary solutions has not been completed yet. At this moment the above ar- 
gument does not seem to apply directly to such 7's. But we believe the method 
developed in this paper can make a contribution to show the instability for those 
7's. We will leave them for future study. 

The paper proceeds as follows. The first half of this article is devoted to develop 
the linear theory: finding the largest growing mode of the linearized Euler-Poisson 
system and deriving some regularity of it. While the linear instability was studied 
in 0, it was done in Lagrangian formulation and it does not give the precise growth 
rate. For our purpose we shall demonstrate the explicit linear instability analysis 
in Eulerian coordinates. In Section 2, we formulate a variational problem to find 
the biggest eigenvalue ^/Jlq and corresponding eigenfunction. In the subsequent 
section, the fast decay property of the largest growing mode is derived. In Section 

4, we show that ^/JIq dominates exponential growth rate of any solutions for the 
linearized system. In the other half the nonlinear analysis is carried out. In Section 

5, the weighted instant energy and total energy are introduced and it will be shown 
that the total energy is bounded by the instant energy containing only temporal 
derivative terms. We perform the weighted energy estimates in Section 6 to get 
the precise estimate of (1.13). Finally in Section 7, the bootstrap argument and 
Theorem 1.1 will be proven. 



Firstly, we study the linearized Euler-Poisson equations for the spherically sym- 
metric case. We assume p > which is our concern. Letting p = p n + a, we 
linearize (1.4) and (1.5) around a given steady state (po,0) and get the linearized 
Euler-Poisson equations in terms of a and u. 



2. Existence of the Largest Growing Mode 



(2.1) 




(2.2) 
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To find a growing mode for (2.1) and (2.2), let a = e xt <p(r) and u = e xt ip{r). 
Then (2.1) becomes 

(2.3) A^+^V + Po^ + -PoV' = 0. 

Or Or r 

Since ^-ip + Poff + ^Pofp = -ps(r 2 Poip)r, (2.3) gives a simple relation between <p 
and ip: 

(2.4) = I <j>{s)s 2 ds. 

r Po Jo 

Similarly, (2.2) becomes 

(2.5) AV, + A^pf 2 ^; + 7(7 - 2)pr 3 |^] + % jf 4>s 2 ds = 0. 
Multiplying (2.5) by and using (2.3), we get 

*± =7Pr 2 (-A0)' + 7(7 - 2K- 3 p[,(-A</>) + r <Ks)da) 

=7Pr 2 boV + 2p ^ + poV/' 



7-2r //„,. , o„/,.' , ..„,." 2 PoV> , 2p'V , 2^0^' ! 



+ 7(7 - VpVp'M + P^' + ^ + ^ 

=ipVr + h 2 P r 2 p'o + %r v + hpr 2 p'o' - 2 ^pV + 2 ^pVp' 

, / 0A 7-3/ /x2 . 27(7-2) 2 , 47rp 0l , 
+ 7(7-2)/^ ( Po ) + pi p + __]^ 

where ' = ^ and this is the 2nd order ordinary differential equation. For the 

further simplification recall that po satisfies (1.6): A-^p'^ 2 p' {) + |f p s 2 ds = 0. 
1 d r 

Ar 2 <2r L 



Compute ^3 ^I 7 " 2 ' (1-6)], and then we get the following relation: 

7-2 // . 27 7 _ 2 , . , n 7-3/ /\2 1 47T/? 

TPo Po + — Po Po + 7(7 - 2)pJ (po) + — ^- = 0. 
In turn we have 

A 2 V> 7-1/// , r 2 7-2 / 27 7-1, ;/ .27(7-2) 2 , 27 - 
(2-6) -T- = 7P0 ^ + [7 Po Po + — Po ]V> + I " Po Po - -2-P0 W- 

J\ T T T 

Multiply (2.6) by ^ll to obtain the following: 

- r 2 p^" + br^rVo + 2r^' + [2(7 - ^p^Po - WW 
(2 ?) = r 2 pli," + (r 2 poOV + 2[2-l r (pl)' - p^ 

= (r 2 p^y + 2[^r(piy-p^ 

Denote the RHS of (2.7) by Lip: 

(2.8) Lip = (r 2 plip')' + 2[^—^r{pl)' p^ 
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Note that the linear operator L is self-adjoint and hence A 2 is real. 

Lemma 2.1. Suppose \ an d u) (lu > Q) satisfy L\ — lu 2 ^-X- Define ip — 
— -^z{ T Pox)' ■ We assume \ an d f are well defined admissible functions in a 
suitable sense which will be clarified later on. Then a = e ±UJt ip and u — e ±UJt x are 
a solution pair of the linearized equations (2.1) and (2.2). 

Proof. This is obvious by the definition of \ an d ¥?■ One can keep track of the 
derivation of L to sec it.D 

Lemma 2.1. tells us that (e^ip, e"'x) satisfying all the assumptions is a growing 
mode for the linearized Eulcr-Poisson equations. Next we show such a growing 
mode actually exists when 7 = §■ This can be done by looking at the eigenvalue 
problem of the operator L due to Lemma 2.1. In other words, we only need to find 
ip and A (A > 0) such that Lip = A 2£ ^-i/>j where Lip is defined in (2.8). 

From now on we fix 7 = | and A = ^ and corresponding 

1 



Po(r) 



(l + r 2 )* 



The starting equation is Lip = /j,^^-ip on (0, 00), where fi = A 2 . It is well known 
that the largest eigenvalue ziq is given by a variational formula: 



po = su P[f^y : W) < < °°] 



where 



Q(iP) = (Lip, iP) = - I r 2 pl(ip' fdr + 2 / 0—^r{pZ)' - P ^ 2 dr 
Jo Jo 7 

/■oo 2 roo o 2 _ 1 



and 



A-f J Q Ay 4tt J (l + r 2 )§ 



Hence once the above formula attains the sup, the largest eigenvalue /z and 
corresponding eigenfunction ip of the linear operator L gives a largest growing 
mode of the linearized equations. In order to carry it out, first define a norm for 

any VeC c °° (0,oo), 

r 2 f°° 2 15 f 00 r 2 

Let H = C£°(0, 00) in the above norm. It is clear that jj^y = j^j for any nonzero 
constant c. Thus the variatonal problem can be rephrased as to find a maximum 
Ho of Q(ip) on H under the normalization condition I(ip) = 1. 



Proposition 2.2. There exists a ipo £ H such that I{ipo) = 1 and Q(tpo) = Mo 
(7*o > 0J i.e. ffce sup fio is attained on H. 
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Proof: First we claim fia > 0. Consider ip = ^Jr. 

Since //o > fjy|) an( i > 0> it 1S enough to show Q{^/r) > 0. 

Q(Vf) _ (- J + 6) ^* - 8 f 
23 r 1 ,„ „ r 1 



4 L 2(l + r 2 ) 2Jo L 3(l + r 2 )3 Jo 
8 3 

q 2 

And note that the positive part of Q is uniformly bounded by I because p^ppys = 
O(i) and — - — g- = 0(A) for sufficiently large r. This implies ^o is finite. To 

(l+r 2 )5 

show fi is attained on H let {ip n } be a maximizing sequence i.e. 

Qi^Pn) / Mo as n — ► 00 and I {ipn) = 1 for all n. 

Let ^0 be its weak limit. Then by the lower semicontinuity of weak convergence, 
we have 

Umillf / (1+^2)4 ^ > / ( 1+ 1 r 2)4 ^. 

„2 /■ 2 

i.2j.n / T .1.2. 



liminf / r ip„dr > / wipt\dr. 

I (l + r 2)!^" -J (l+ r 2 )5 ° 



Claim 1. (Compactness of the positive part) There exists a subsequence {ipn k } 
of {ip n } such that 

f r 2 f r 2 

J ( 1 + r 2 )4 < rfr — » y (l + r 2)4 ^0<*r aS "* — °°- 

Claim 2. /(^o ) = 1- 

Since Claim 1 immediately implies Q{ipo) — Mo, the conclusion follows from 
Claim 1 and Claim 2. It remains to prove Claim 1 and Claim 2. Claim 2 follows 
from a simple scaling argument. Suppose |f / - — r 2 ^s ip 2 dr — a 2 < 1. Then 

i( — ) = 1 and Q{ — ) = —fi Q > ^ 

a a a z 

which is a contradiction to the definition of /io- To prove Claim 1, first observe that 

roo 2 1 r°° ^2 



r°° r 2 1 roc r 2 

Jr JTTW^- (l + R 2 )d R 0T^ €dr ~ 15(l + i? 2 )§ 
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Fix R > 0. On the finite interval (0,i?), since Tj-ppiyc^dr ~ llV'l^^^g)), 

^ (iS)^ 2dT ~ I^Hi 2 (s R (0))' and /cf TTWW'') 2 ^ ~ IIV ,, |||2 (Br(0)) where 
B R (0) is a ball with radius i? in R 3 , we can apply the Rellich-Kondrachov Com- 
pactness theorem, which says H 1 (Bp(0)) is compactly embedded in L q (Bn(0)) for 
each 1 < q < 6, . So there exists a subsequence {Vw} sucn that 



o (l + r 2 )^ n " J a (1 + r 2 ^ 



« r 2 



— > / — ^rr%dr, as n fc — ► oo. 



For any given e > 0, choose i? > large enough so that - — 3- < e. 

15(1 + -R 2 )2 

Now take the limit of — ► 00. Since e > is arbitrary, this finishes the proof of 
Proposition 2.1.D 

Finally, let us make sure tpo is in fact an eigenfunction corresponding to hq. 
Consider a perturbation ipo + er\ around t/> . Since fi = Q(ipo), 

-J7-— — r < Q{i>o) = Mo 
I(ipo + en) 

for all sufficiently small e and all admissible function 77. And hence we have 
Q(ipo + eq) < + eq) 

2 

QWo) + MLtPo, V) + e 2 Q(v) < Q(ih){I(ih) + M—iPo, v) + e 2 I(v)} 

7 

Here (-,-) is the standard inner product in L 2 (0, 00). Note I(ipo) — 1 an d Q{v) < 

fioI(ri). Hence in order for the above inequality to hold for all e and 77, the coefficient 

2 2 
of e should vanish, i.e. (LifroiV) = ^oi^Tf^OiV) f° r all 77. Thus, Lip a = /J-o^rr^o 

and therefore y/JIo and ipo give a largest growing mode for the linearized Euler- 

Poisson equations. 



3. The Regularity of the Largest Growing Mode 
Recall tpo satisfies the following 2nd linear ordinary differential equation: 

!5mo r 2 . r 2 „ 3r 2 - 1 , 

— (i7^° = ( (TT^ o) + 2 (TT^° 

(3.1) r 2 ,„ 2r-4r 3 ,, „ 3r 2 - 1 , 



The interior regularity easily follows from the elliptic theory of the 2nd order 
differential equations: since each coefficient of ipQ, tp' and tpo in (3-1) is in C°°(e, R) 
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where e is a small enough positive fixed number and R is a large enough fixed 
number, ipo is also C°° on (e,R). In the following two subsections we investigate 
the behavior of i/'o when r is either very small or very large. 



3.1. The behavior of ?/>o near the origin. In this first subsection we show tpo 
is analytic near the origin. This rather surprising property easily follows from the 
classical theorem by Frobenius from the ODE theory. For reference, see ■ Before 
we prove the analyticity, it will be shown that the maximizing property implies the 
boundedness of ipo at the origin. 



Lemma 3.1. There exists a decreasing sequence {e^} with \ such that 
V'o(efc)V , o(efc) > for each k. 

Proof. Consider tpi — Qipo, where O is a lipschitz cutoff function defined by 

'0 ifO<r<e, 
6(r) = { r-e if e < r < 2e, 
e if r > 2e. 

Here, e is a sufficiently small positive number to be clarified. Then clearly tpi is in 
the admissible set H. First note that I(ipi) < I(eipo) since ipf < e 2 ipQ. Next let us 
look at the difference between Q(ipi) and Q(eipo). Since ipi has the same value as eipo 
onr > 2e and ipf = ipo^(e,2e) +2(r- e)V'oV'o 1 ( e ,2e) + (r- e) 2 ^ 2 l( £ ,2 e ) +e 2 ^o 2l (2 e ,oo), 
Q(V>i) - Q(eV'o) is 

(3.2) 

' (l + r a )3 [^ + ( r - £ ) 2?/, o 2 + 2(r - e)M' ]dr + 2^ ( f + ~ 2 | 4 ( r - e )^o* 

2e r 2 /-2£ o 2 _ i 



/o (1 + ^ 2 ) 3 " U ^0 (1+^ 2 ) 4 

After collecting similar terms, (3.2) can be written as following: 

(3.3) 

f 2e r 2 r 2e 2r 2 (r - e) 

repair + r (2e-r)^ 2 dr- / J ^dr 



Q (l + r 2)3 ™ j ( (i + r 2)3 v _/ £ (l + r 2)3 

2 - 6r 2 2 2 l , [ 2t 5r 4 - 12er 3 - 3r 2 + 4re 



o (1 + ^ 2 ) 4 ™ J* (l + r 2 ) 

Call the first three terms of (3.3) (I) and the last two terms (II). Now suppose the 
lemma were false. Then for any small enough e > 0, we may assume i/'oV'o < 011 
(0, 2e). It is obvious to see that (I) is positive. And because ipQ is decreasing on 
(0,2c), 

2-6r 2 f 2e 2-6r 2 ,,2, 



(1 + H)4 ™ y £ (l + r 2) 
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and hence 



2£ 2 - 6r 2 „ , , , f 2e 5r 4 - 12er 3 - 3r 2 + Are 



2£ 5r 4 - 12er 3 - 6e 2 r 2 + e 2 , f 2e e 2 + Aer - 3r 2 2 

/ — t: — %dr 



(1 + r 2 ) 4 .,, , . , , , 

= (III) + (IV). 

(Ill) > since 5r 4 - 12er 3 - 6e 2 r 2 + e 2 > on (e, 2e) for sufficiently small e > 0. 
Claim. (IV) > 0. 

Claim implies Q(ipi) - Q(etpo) =(I)+(II)> 0. This leads to 

Qbh) Qi&M = QW = 

/(^) > /(e^o) /(V'o) 
which is a contradiction to the definition of fiQ. Therefore, to finish the proof it 
suffices to verify Claim. Let a be a zero of e 2 + 4er — 3r 2 = with e < a < 2e. 
Then 

e 2 + 4 e r_3r 2 = ^° ife ^<a, 
[<0 ifa<r<2e. 

The decreasing assumption on ip 2 will be again used: 

ib 2 (r) f a ib 2 (r) f 2e 

(IV) > min , ytU ' A / e 2 + 4er - 3r 2 dr + max , „\, / e 2 + 4er - 3r 2 dr 
(1 + r 2 ) 4 J e [ a ,2e] (1 + r 2 ) 4 J Q 



V'o(Q') 
(1 + a 2 ) 4 

0. 



,2e 

/ e 2 + 4er - 3r 2 dr 



The first equality holds because is decreasing and the second one is simply 

due to the fact J £ e 2 + 4er — 3r 2 dr — 0. This completes the argument. □ 

The Frobcnius theorem with Lemma 3.1 gives rise to the following analytic prop- 
erty of ipo around the origin. 

Lemma 3.2. ipQ is analytic at r — and moreover ipo(f) = ar + o(r 2 ) around the 
origin where a is a constant. 

Proof. In order to employ the Frobenius theorem, we need to show the equation 
(3.1) has a regular singular point at r = 0. To check this out in the context of |2] 
(p. 215), we rewrite (3.1) in the following form 

,// 2-4r 2 ,, r 23r 2 -l 15^0,., 2^1/ 
r(l + r z ) r z 1 + r z An 

Let P(r) and Q(r) be coefficients of ip'o and V'o respectively in the above. Then 
it is clear that rP(r) and r 2 Q(r) are analytic at r = 0, which means r = is a 
regular singular point of (3.1). Let po, g be the zeroth order term of rP(r) and 
r 2 Q(r) respectively. It is easy to check pq — 2 and qo = —2. The indicial equation 
r(r — 1) + por + qo = has two roots r\ = 1 or r<i = — 2. Hence by the Frobenius 
theorem ipo has a power series representation of either yi(r) = r^*L Q a n r n or 
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V2{r) = at/i (r) In r + r 2 X)^Lo^ nr ™- However, 2/2(7") is impossible by Lemma 3.1 
and therefore tpo — rJ2^=o a nr n is obtained. □ 



3.2. Asymptotic behavior of ipo. In this subsection we first observe V'o( r ) i s 
nonincreasing near the infinity owing to the maximizing property and then the fast 
decay of tpo in an appropriate sense will be shown by using a standard bootstrap 
argument. 



Lemma 3.3. There exists an increasing sequence {Rk} with Rk /* 00 such that 
i>o{RkWo{Rk) < for each k. 



Proof. Consider tpi — Qtpo, where 6 is a lipschitz cutoff function defined by 
0(r) = < 



1 if < r < R, 

R+l-r if R<r < R+l, 
if r > R+ 1. 



Here, R is a sufficiently large number to be determined. Note that (ip[) 2 = 

W)) 2 l(0,fl) + («+!- ^) 2 (^ ) 2 l(K,K+l) + V>ol(fl,H+l) - 2(^+1-0^01(^+1)- 

Then clearly Vi € H. Let us compute /(^i) and Q(ipi). Recall I(ipo) = 1 and 
Q(V'o) = Mo- 



1C /-OO 2 ir /-fl+1 



47T 

il-Cfl 



[/ (1 

JR. 



(1 + r 2 ) 
(i?+l-r) 2 



(l + r 2 )§ 



-■00^ - 



Q(Vi) =Q(Vo) - [ 

r--R+l 



(l + r 2 )§ Jr+i (1 + r 2 )" 

3r 2 - 1 



„ (TT^^ + 2 A (T+^yi^ 



/« (1 + r 2 ) 3 

[R+l 

-2 (R+l-r) 
Jr 



^alr + ^ (R + l- rf (l + r2) ^ dr 



R+l 



(TT^T^ dr - 2 J R {R + 1 - r) (1 + r 2 )* 



3r 2 - 1 



11 ' : r 2 



>Mo-[-2/ (.R+l-r) 

! / (l-(i? + S -r) 2 ) * 

JR 



R+l 



r 



+ 2 



„ (1+r 2 ) 
00 3r 2 - 1 



3^0* 



+ 2 



o 2 dr] 



/■fl+i 

=Mo + 2 / (R+l-r) 

JR 



(1 + r 2 ) 3 



ipoi>'o dr - d r 



Now suppose the proposition were false. Then there exists large enough Rq such 
that if R > R , then VoV'o > on (R, R + 1). On the other hand, D R < §C R for 
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R > Rq > 0, where K is a fixed constant. So we get the following inequalities. 

QOi) Mo + 2 J^ +1 (R + 1 - r) (1+ r r 2 a)j VoV'cA ~ D R 
~ 1-Cr 
^ Mo - \Cr (hq - ^)C R 

> ^Cr~-^ + 1-Cr 

Choosing sufficiently large R > Rq > so that ^ > we § et Tj^jrfj > Mo which 
contradicts the definition of /Uo- Therefore the lemma follows. □ 

Lemma 3.3 with the integration by parts leads to the next proposition that shows 
the fast decay of ipo in a suitable sense. It will play a key role when we prove the 
main theorem in Section 7 in a sense that the initial data with different weights are 
not essentially different. 



Proposition 3.4. J °° ^ + ^ iffidr and J °° ^p^js ^odr are bounded for each non- 
negative integer n. 

Proof. Let Rk be given as in Lemma 3.3. Multiply (3.1) by r n ip and integrate 
over (0, Rk). Then we get 
(3.4) 

r-Rfc r n+2 



15 M o f nk r n +-' . [ R * r 2 f R " n 3r 2 - 1 2 

The first term on the RHS of (3.4) can be rewritten as following using the 
integration by parts, 

r n+2 rRk ~2 

Mo)o k - / (nr n -^o + r n ^) 7 —-^^dr 



l( 1 + r 2)3^0^J0 ^ V- ru , ■ ru/ (1 + r2)3 



fn k n+1 

f>n+2 E n+l r Rk r n+2 



fl* ( n _ 5)r™+ 2 + (n + l)r 

2j Q {iT^V 



ipfidr 



Plugging this into (3.4), we get the following relation. 



Rk r n+2 fRk n+2 

R n+1 n 

(3-5) + {l + Rlf [ 2^ {Rk) R ^'o(Rk)MR k )} 



n 



2T3 — ^ + 2 / Vo 2 ^ 



2 7o (1 + ^ 2 ) 4 rU Jo (1 + ^ 2 ) 

Observe that for each Rk, each term of the LHS of (3.5) is nonncgative be- 
cause of the previous lemma and the small r parts of the RHS of (3.5) are finite 
due to the behavior of tpo near the origin. First, when n = 0,1,2,3, all terms 
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in the RHS of (3.5) are finite for any Rk since they are uniformly bounded by 
if Jo°° — ~ — z~^od r = 1. As taking the limit of Rk , we know each term of the LHS 

converges as Rk — ► oo. In particular, J"°° — rip^dr and f°° -rrnpzw (V'o) 2 dr are 

bounded for n = 0, 1, 2, 3. The standard induction on n with (3.4) gives the desired 
result for general n.D 



4. The Linear Growth Rate 

In this section, we show that ^/po is the dominating exponential growth for the 
linearized Euler-Poisson equations. The linearized Euler-Poisson equations (2.1) 
and (2.2) for 7 = | are 

(4.1) a t + -^(r 2 Po u) r = 



47T , -$ N 47T 

15 



(4.2) Ut + —(p Q 5 a) r + — I as 2 ds = 



Multiply (4.2) by ||r po, take the t derivative and use (4.1) to get rid of a. 

^r 2 p u tt ^-r 2 p (p^(T t ) r -15p J a t s 2 ds 
1 -i 

=r Po(~^P 5 (r 2 pou)r)r + Vor 2 p 2 Q u 

J3 JL 1 

=r 2 p^u rr + {r 2 p (p^ ) r + p$ (r 2 p )r}u r 
1 

-2 



+ {r 2 Pa(— Po 5 (r 2 po)r)r + ^r 2 pl}u 



After putting back p — — 5- in the above, we get an equivalent 2nd order 

(l+r 2 )3 

equation for u = ^ : 

15 r 2 r 2 3r 2 — 1 

(4 - 3) = ( (TT^F^ + ^(TT^* 

Here we use ^ instead of u in order to distinguish the linear analysis from the 
nonlinear one. In this way of writing it is also easy to compare (4.3) with (2.8). 
Note that \E" in (4.3) is a function of both t and r while tp in (2.8) is a function of 
only r. Denote the RHS of (4.3) by Lo^ ■ L is basically same as L in Section 2. 
Define the following quantities. 

™ r ^ 15 2 15 
VVo(r) = -r-r po 



47T 47T (l+ r 2)2 

WfWiv = (Wf, f) = / Wf 2 dr where W is a given weight function 
Jo 

r°° r 2 Z" 00 1 
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Lemma 4.1. For every solution ^ to (4-3) there exists Cp, a ,C^ a a > such that 
(1) 



ll*(*)lko , ll*t(*)lk„ < C^e^(\\^ t (0)\\ Wo + \\*(0)\\ Wo + y/l%(p)). 
(2) For any a > 1, 

||af +1 *(t)lk < ^ 0;ae v^*(||^(o)|| Wo + ||*(o)|k + v^mo)) 



+c AJO , Q ^(||aj +1 *(o)|| Wo + yJp ai *(o))- 

i=l 

Proof. Take the inner product of (4.3) with Then we get 

(Wo *«,*«) = (£o*,*t) 
^^(W *t,*t) = ^(io*,*). 

The above equivalence comes from the self-adjointness of L . Next integrate the 
above with respect to t to get 

(4.4) (Wb* t (t), * t (t)) = (L *(*), *(*)) + (W * t (0), * t (0)) - (L o *(0), *(0)). 

Since (L *(t), *(*)) < /xo(W tf (t), *(t)) for all * and -(L tf (0), tf (0)) < P*(0), 
from (4.4) we get 

(4.5) \\Mt)\\ 2 w < ^(t)\\ 2 Wo + \\*t(0)\\ 2 w + *V0). 

Since < J * 11**0)11^^ + ll*(°)llw> Pegging this into (4.5), we get 

||*t(*)lko <^f ll^tWHwbdr + C(||* t (0)|k o + ||*(0)|| Wo + v^MO))- 
Jo 

By Gronwall's inequality, we obtain 



|*t(*)||wb < Ce^dl^^O)!!^ + ||*(0)|| w „ + y/P*W) and in success 



||*(*)lko < Ce^iWMVWwo + ll*(0)lko + v^MO))- 

Notice that C only depends on fj, . For higher derivatives, take df of (4.3): 
W (r)d? = L d?V. Take the inner product of this with t d t a * to get 

R Q+1 *(*)llk 

= (w a t Q+1 *(t),a t Q+1 *(i)) 

= (LoW^W*)) + (T^ 9r +1 *t(o),a« +1 * t (o)) - (i a t a *(o),a t Q *(o)) 

< f, \\d^(t)\\ 2 Wo + \\d? +1 *(0)\\w o + Pe?*(°)- 
Thus (2) easily follows. □ 
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Next we show that the energy estimates with Lemma 4.1 lead to the same ex- 
ponential growth rate on the a satisfying (4.1) and (4.2). To avoid the confusion 
with the nonlinear analysis, we use $ instead of a. 

Define a weight function for $ by 

V (r) = r 2 p^ =r 2 (l + r 2 ) 2 . 

Notice that (Vo,Wo) is chosen not randomly but to be a symmetrizer of (4.1) 
and (4.2) that makes the energy estimate work. ||$||vb + ||^|| w resembles the real 
energy (1.7). 

Lemma 4.2. There are constants C\,Ci > such that 

\\d?m\\v < c^dr^mwo + wm\w + h wo)n Wo) j or an ya >i 

and \\<f>(t)\\ Vo < C 2 eVP5-*(|| **(0) ||^ + 1|*(0)|| vb + II *t(0) || w + 1| *(0) || w + V^(0))- 
Proof. Fix a > 1. Compute the following equation: 

r°° 4 1*5 
J [p • S«(4.1) + ^Pod?* ■ d?(4.2)]r 2 dr = 0. 

_ 4 

The choice of weight functions p 5 and jfpo above yields a nice cancellation after 

integrating by parts, i.e. pg*d?$ ■ {r 2 p d?^) r + r 2 p a d^f> ■ (p ^ d?<Z>) r dr = 
and it results in 

(4.6) 

2Jt I bo*(9t*) 2 + ^Po(d?*f]r 2 dr + 15 J p d?*( j d?<f>s 2 ds)dr = 0. 

On the other hand, by (4.1), we get J Q r d?<S>s 2 ds = -r 2 po^"" 1 * and hence (4.6) 
can be rewritten as 
(4-7) 

As taking J* of (4.7), the desired result is obtained. Note that we have used 
pl < Po- As for a = 0, utilize ||$(t)|| Vo < /* ||$ t (r)|| Vo rfr + ||$(0)||y o .D 

The next two lemmas show that /zo also determines the exponential growth rate 
even with strong weights. Lemma 4.3 will play a crucial role in the proof of The- 
orem 1.1. Main idea of proofs is to utilize the linear operator L. The results only 
contain if? estimates and the estimates on $ can be derived similarly. 

Lemma 4.3. For any a > 0, there exists C^ > such that 

poo 2 poo -i 

where I a = ^7=0 ll^t^Wllwb Sf=!o ^9|*(0) * s given initial data, for every 
solution * to (4.3). 
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Proof. Multiply (4.3) by df ^ and integrate to get 
(4.8) 

i-oo roc 2 /-oo o 2 _ i 

If we move two negative terms in the RHS of (4.8) into the LHS, we obtain 

roc 2 roo i 

foe „2 roo 

= 6 / n I 2^4 {d^fdr- / W^+^vMr 
Jo U + r J Jo 

< 6 I J^yJ^dr + 2(ll^ +2 *Wll^o + II W*)llvr ) 

<c(\\a? + H(t)\\* Wo + \\d?*m 2 w ) 

But the last quantity is bounded by Ce 2 ^' t L a by Lemma 4.1 and the conclusion 
follows. □ 

For the next lemma define the weight functions Wi for I < by 
Wi{r) = (1 + r 2 )-^W (r) = -^r 2 . 

47T 

We remark that Wi is a linear version of symmetrizers Si defined in Section 5. 

Lemma 4.4. For any a > 0, there exists C Mo > such that 

\\d?*{t)\\w k < C^e^y/i;, 
where I a is defined in Lemma 4-3, for every solution W to (4-3). 

Proof. Let us only consider a = 0,1. Other cases can be treated in the similar 
way. Letting L k = (1 + r 2 )?L , wc get (W fe = (L fe *,*t). But because 
Lfe is not self-adjoint any more, we do not have a simple equivalent expression as 
in Lemma 4.1. Instead, we try to find a relation between Lk and L$. The follow- 
ing identities are needed and they are obtained by definitions and straightforward 
computations. 

(*(l + r 2 )4) r = (l + r 2 )H r + ^r(l + r 2 )*- 1 * 

r°° i. r 2 f°° i 'V 2 — 1 

(L fc *,*) = jf ( l + r 2 )^(___v &r)r(ir + 2y o (i + r^—-*"* 

/>oo 2 

(L (*(l + r 2 )*),*(l + r 2 )*) = - / {(1 + r 2 )H 2 + fcr(l + r 2 )*" 1 *^ 

+ ^r 2 (l + r 2 )|- 2 * 2 }dr + 2 f° ^-^ (1 + r 2 )t * 2 dr 
4 7 (l + r 2 ) 4 ' 

The last two identities imply the following: 

(i fc *,*) = (L (vI/(l + r 2 )*),vl/(l + r 2 )4) + - jf ^-^(*(1 + r 2 )4 ) 2 rfr . 
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Define Qk and Ik similar to Q and I in Section 2 as following 

r°° r 4 

Q fc (*) ee (L (*(l + r 2 )*), *(1 + r 2 )^ ) - /J ^ (*(1 + r 2 )*) 2 dr, 

!*(*) = /(^(l + r 2 )*^^*,*). 

where /3 is a small positive constant. Then by doing the same variational analysis 
in Section 2, one can show that there exists fj,k > such that /ik is the maximum 
of a functional j^t§j and hence Qk(^) < Hklk{^)- 

Claim, [ik < Mo- 

To see the Claim, pick a <f"i such that Qk(^i) = and Ik(^i) = 1- By the 
definition of Qfc, 

Q fc (*i) = (L„(*i(l + r 2 )*), ^(1 + r 2 )^) - /? / ———(^(l + r^fdr. 

Jo (i- + r y 

Since (i (*i(l + r 2 )4),*i(l+r 2 )4) < Mo /(*i(l+r 2 )*) -Mo4(*i) =Mo, 



Thus, 



(4.9) 



Ofc(*i) - Mfc < Mo - /? / T 2^5 + ^ 2 )~) 2 ^ < Mo- 
Jo {i + r )' 



1,2 rco 4 

(L fc *, *) = QkW + (/?+—) y o (1^5)6 (*a + r 2 )* 

L2 /-ex) 4 
< *) + 03 + -) / -ToTF* 2 * 

4 Jo fl + r 2 )^- 



/o (1 + r 2 ) — 
Now we are ready to go back to {W k ^>tt, *t) = (£fc*, *t)- 

„2 



/■oo r 2 /-oo r 3 /-oo 3? ,2 _ j 

= / 7i -^r^rVrtdr-k _ 1 - 5 -M t dr + 2 / 

Jo (1 + r 2 ) — Jo (1 + r 2 )— Jo (1 + r 2 ) — 

1 d f 00 /*°° ^ 

= 9 3l( L fe*' *) - * / (- -^-)r**tdr - fc / ^(fr- 
iar Jo (1 + r 2 ) — J (1 + r 2 ) — 

Id r°° r 3 



Hence, 

l_d 

OO 



(4.10) 



k / : —^r^tdr. 

Jo (l + r 2 )V 



Let k = 1. By Lemma 4.1 and 4.3, using the Cauchy-Schwartz inequality, 

/■oo 3 poo 3 

~ / ( n ^ 2 , T )r^tdr - / r 7 *MMr < Ce 2 ^I 

Jo (1 + r 2 ) 2 Jo (1 + r 2 ) 2 

for a constant C and initial data I . Rewriting the above (4.10) when k = 1, get 
\^W^ t )<\^^ ) + Ce^I . 
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Taking the integral with respect to t and using (4.9), we have 
(Wi* t , * t ) < (£1*, *) + Ce 2 ^ 4 ^ 




1 Z" 00 r 4 



* 2 dr + Ce 2 ^*7 ( 



o 



(4.11) (Wi*t,*t) < /ii(Wi*,*) + CeV«*Jo. 



Since < / * ||* t (r)|| w dr + ||*(0)|| w , combining this with (4.11), get 



Jo 

Since we have 1//I1 < y^to, Gronwall inequality gives H^tJHwi < Ce^^^fT^ as 
well as ||*t(<)||wi < Ce^ l \fh. The standard induction on fe with (4.10) claims 
the desired result for all k.O 



5. Weighted Instant Energy and Weighted Total Energy 

In this section, by the utilization of symmetrizers of the Euler-Poisoon system, 
we introduce suitable measurements of perturbations a and u of steady states re- 
sembling weighted Sobolev norms: the weighted instant energy Ei and the weighted 
total energy where I is an index associated to weights. The symmetrizers will 
play the same role as weights in the linear analysis for small solutions. Then the 
total energy is shown to be bounded by the instant energy under a certain smallness 
assumption in using the equations directly, which makes it sufficient to play only 
with the instant energy. The weighted energy estimates for the instant energy will 
be carried out in the next section. Before going any farther we remark that it is con- 
venient to work on rectangular coordinates rather than polar coordinates because 
one can avoid the singularity of the origin coming from the spherical symmetry. 

We are interested in sufficiently small solutions a, u satisfying the neutrality 
condition 



For such small solutions the Euler-Poisson system (1.1), (1.2) and (1.3) can be 
rewritten in the rectangular coordinates as the following: 



47T 4 107T —2 

(5.3) u t + (u • V)u + — {p a + «7)~*V<t - {^r-p 5 v + h{a,p )}Vp + W = 





(5.2) 



0* + (Po + 0-)V • u + V(p + <j) ■ u = 



(5.4) 



A<f> = 4na 
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where h(a,po) = — ^§{(/Qo + o") = — Po 5 + f Po '' a } represents higher order terms. 



u takes the vector form, i.e. u(x, t) = u{r, t)f, where x G I 
each component of ti by Note V x u = 0. 

Now let us consider the symmetrizers Si where I e 
system. 



Si 





(po+a) 1+; 








(po + <?) 1+l 




I a; and we denote 



for the Eulcr-Poisson 



(po + ^) 1+ 7 



Define the instant energy £i (t) and the total energy £i (t) by 

3 



£,(t) = V / diu\diu 2 , ■ (dia, dlu\ dlu\ %u 3 )dx 

j=0 
3 

3 

^w^EE / ^(Po+a)-^^iar^i 2 + (Po+a) i+ '+t| 9r ^ U |^ 



J = i=0 

j=0 i=0 

Here 9 X represents any spatial first derivatives. Note that ^(t) is a part of the 
real energy 2E of (p, u) defined in (1.7). The case I = 0, however, is not enough for 
our purpose because we cannot close the energy estimate at the step I = 0. Being 
convinced by Lemma 4.4, we try different Vs as well. As I < gets smaller, the 
weights become stronger due to the behavior of po- Unfortunately, as I varies, new 
quadratic terms come out while performing the energy estimates. This phenomenon 
seems undesirable but it turns out that they are equipped with weaker weights. And 
that opens another door. 

Observe that the weights of mixed derivative terms in Si are different, in fact a 
little better, from the ones of temporal derivative terms. This is not a coincidence 
but rather a nature of the system; the same feature can be seen in the linear 
analysis. See Lemma 4.3. £i contains all the spatial and mixed derivatives of a, u 
and it is easy to see that 

3 
3=0 

Now we want to show the converse, in other words, £\ is also bounded by £i under 
a certain smallness assumption: 



(5.5) 



£_b +£_i +£_s < 

5 5 5 
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where 9\ is a sufficiently small constant. In order to appreciate the utilization of 
such an assumption, first we prove the next lemma. 

Notation. J dx represents J R3 dx and when dealing with line integrals / dr, each 
end value will be specified. 

Lemma 5.1. Suppose (5.5) holds for < t < T. Then there exists a constant 

C > such that for each < t < T, 

(5.6) 

sup l^i + l^bl + i ? + "vl ' + I , !V | + |Vm| ^ 

iei 3 Po + cr po + cr (po + cr) io (p + a) io (p + a) w 
In particular, the assumption (5.1) is justified. 

Proof. The above smallness assumption (5.5) together with the Sobolev imbed- 
ding theorem yields the result. To see how it works, let us apply the Sobolev 
imbedding theorem to I ff I. 

° 1 po+cr' 

(5.7) 



S up|^| 2 <C||^||^ (R , 

Po + cr 7 po 



< C { J \-^—\ 2 dx + I \-^\ 2 dx + / l-^-l 2 ^' 



+ [\ ° 9x<yPa+(T \ 2 dx I [\ dx<1 dx ^ Po + <7 " > \ 2 dx 

J P0 + cr po + cr J P0 + O- po + cr 

| ° d * ipa+a) fdx+ /{l-g-H^ + ^l^cfa} 

P0 + cr po + cr 7 po + cr p + cr 

Recall the behavior of p , namely po(r) = 0(r~ 5 ) for large r. Hence | | and 

d 2 

| | arc uniformly bounded. Thus (5.7) becomes 

cr l0 r< o.o 1.1 2,2 % 

sup 2 <C{(£_6 +£_r +£_ s 

po + cr 5 5 5 

CT l0 , 0,0 1.1 2.2 s 

+ SU P (£_£ + £-1 + £-2 

po + cr 5 5 5 

( 5 ' 8 ) , , Cr 2 ^(T 2 —1,1 

+ SU P | — sup |- -j- £_7 } 

P0 + cr (p + cr)To 5 

<C(0 1+ su P |^— l^+supl-^-^supl— ^-t-I^i). 

P0 + cr Po + cr (p + cr) To 

We can get similar estimates to (5.8) for other terms in (5.6). Call the LHS of 
(5.6) S. Consequently, we obtain the following: 

(5.9) s 2 < ce l + cs 2 e 1 + cs i e l 

Since 9\ is small enough, (5.9) immediately implies the lemma.D 

Since £i = tf , we only need to show that t\ where 1 < i < j < 3 is bounded 
by C£i; the precise statement is given in the following lemma. Main idea for ac- 
complishing the goal is to estimate the spatial and mixed derivative terms directly 
from the equations in terms of temporal derivative terms. 
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Lemma 5.2. Suppose (5.5) holds for < t < T. Let I < 0. Then there exists 
C > such that for each 1 < i < j < 3 and for < t < T, 

ir(t)<cj2^(t) + cJ2 £ i +i - 

fe=l fe=0 

Proof. Let us start with one spatial derivative terms corresponding to i = 1. 
No temporal derivative terms i.e. when i = 1, j = 1 are treated carefully since it 
is instructive and other cases can be easily shown from it. Various case of j for a 
fixed i will be done in turn. Then we move onto other cases of i. First of all, solve 
(5.2) and (5.3) for Vcr and V • u to get 

15 , „ 4 r 167T -2 „ , 

(5.10) Va = - — (pa + <j)-{u t - —p °Vp cr + V<j>+(u- V)u - hVp } 

47T 75 

(5.11) V-tt= — {o-t + Vpo -it + Vcr-u}. 

Po + o- 

Notice that the estimate on V • u is enough for d x u since Vxu = 0. In order to 

— 1,1 3 1 

get a right weight of £i for a and u, multiply (5.10) and (5.11) by (p + er)~io + 2 

and (po + cr)t+3 respectively, square and integrate them over M 3 . 

(5.12) 

J ( Pa +cr)-i+ l \Va\ 2 dx + J ( Po + *)t +i |V • w| 2 ^ 

< C j(p + cr) 1+ '{| Ut | 2 + \p~ § V P0 | 2 T 2 + |V^| 2 + |(« ■ V)u\ 2 + |Vpo|^ 2 }^ 
+C j ( Po + *)-i +l {\* t \ 2 + |Vpo • u\ 2 + \Va ■ u\ 2 }dx 

<C j( Pa + o-) 1+l \u t \ 2 + ( Pa + a)-i+ l \a t \ 2 dx 

+C J (po + <j) l+l \p hp \ 2 o- 2 + (p + <7)-t+<|Vp • u\ 2 dx 
+C J (p + a) 1+l \(u- V)u\ 2 + (p + a)-^+ l \Va-u\ 2 dx 

+C J (po + a) 1+l \Wp Q \h 2 dx + C j(p + a) 1+l \V4>\ 2 dx 

We rearranged terms according to the order of algebraic degree. The key difficulty 

lies in the potential part and so it will be done last. The first integral is bounded by 

_ 9 3 
C£\ by the definition. Since \p 6 Vpol = 5r(l + r 2 ) < C(p + o-)~~s and similarly 

|Vp | <C( Po + o-)%, 

J (p + o-) 1+l \p~' i \7po\ 2 o- 2 dx + j (po + o-)-7- +l \\7p -u\ 2 dx 
( 5 - 13 ) <cJ( Po + o-)-r,+ l o- 2 dx + C J(po + <y)i+ l \u\ 2 dx 
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For higher order terms, we use the smallness assumption to get the desired 
estimates: Lemma 5.1 is used. 

J (po + <yf +l \{u ■ VH 2 + (po + <j)-r,+ l \Va ■ u\ 2 dx 
(5.14) < / |- ^-^| 2 {(po + <J) I+1 \V ■ u\ 2 + (po + a)-r.+ l \V<r\ 2 }dx 

J {po + <y) 10 

< C0! [ (po + a)f+'|V • u\ 2 + (po + a)-i+ l \Va\ 2 dx 



Hence the higher order terms in the third integral of the RHS in (5.12) can be 
absorbed into its LHS, since 9\ is sufficiently small. Because h is a higher order 
term depending on only po and cr, after applying Lemma 5.1 to the fourth integral 
in (5.2), we get 

(5.15) J ( Po + a) 1+l \V Po \h 2 dx < C#i£ ; ° + 3. 

In order to handle the potential part, we write the poisson equation (5.4) for the 
spherically symmetric case in polar coordinates as 



+ -<p r = iira where 4> r = — ^ / <rs 2 ds. 



}2. 



o 



Recall the L p estimates \\d 4>\\l 2 (r 3 ) < C||c||l 2 (k 3 ) an d in the spherically symmetric 
case it implies that since Y^i,j=i(®idj ( t > ) 2 = 0rr + 

||-0r|U 2 (M3) < C||(T|| L 2( R 3). 

We consider two cases: I > — | and I < — |. In the first case, (p + a) 1+l r 2 ~ 

(1 + r 2 )~ ( i ] r 2 is uniformly bounded, and therefore we get 

(5-16) 

J \ Po + <j) 1+l \\70\ 2 dx = 4tt J~ \ Po + a) 1+l ^ydr < C J~ \^-fr 2 dr < C J a 2 dx. 

The L p estimate has been used at the last inequality. When I < — |, we divide the 
integral into two parts. 

(5-17) 

J (po + <j) 1+l \V4>\ 2 dx = 4tt jf°(po + <J) 1+I tfr 2 dr 

= 4nJ (po + <r) 1+l (f> 2 r r 2 dr + 4ir J (p + <f> 2 r r 2 dr 
= (I) + (II) 

In the unit ball we can do the same as in (5.16), since the weight function is 
bounded. 

(5.18) {I)<C J (y) 2 r 2 dr < C a 2 r 2 dr < C J "(p + a) m a 2 dx, for any m. 



NONLINEAR INSTABILITY IN GRAVITATIONAL EULER-POISSON SYSTEM FOR 7 = § 23 

For the second term (77), we use the neutrality condition J R3 adx = which is 
equivalent to J Q r <rs 2 ds = — <rs 2 ds in polar coordinates. 

(5-19) 

(77) <cj^(l + r 2 )-^ (1 £ as 2 ds) 2 r 2 dr 

-^-^{il + r 2 )-^ J r „ S 2 ds} 2 dr 

POO i rOC 



/•oo i /■ oc roc 2 

pOO 1 P OO 2 /-oo 

<C j ( Pa + cr)-i +l a 2 dx 
From (5.16), (5.18) and (5.19) we conclude that for any I < 0, 
(5.20) J (p + cr) 1+l \Vcj)\ 2 dx < C£° + s. 

Thus, from (5.12), (5.13), (5.14), (5.15) and (5.20), we obtain 



Next we focus on one spatial, one temporal derivative terms, namely the case 
i = 1, j = 2. Take d t of (5.10) and (5.11): 

(5.21) 

_ 3 (Jt , ,4 167T —2 _ , , , _ , 

\<J t = (po + cr) 5 {u t - —— p 6 Vp o- + V0 + (ti • V)u - Wpo} 

7T po + °" '5 

15, .4 lQlT -2 

- i~(Po + 0") 5 W - -^rPo ° Vpo^t + V0 t + (u t • V)u+ (u • V)u t - fttVpo} 

47T 10 

V • tit = — — — {c* + Vp • U + V(T • u} 

(5.22) Po+<7p + <7 

{cr 4t + Vp ' Ut + V(T t ■ U + V<7 • M t }. 

PO + O" 

The first part of (5.21) and (5.22) is bounded by |-2£_Vct| and |-^V • u| that 
have been already estimated at the previous step. Note that | °^_ a | is small enough 
and it does not cause any trouble. The other part has the same structure as before 
and therefore we can do the same: multiply the same weights used in the previous 
case, square and integrate. The potential term is easily taken care of, since the 
dynamics of V<fi t gets simpler and better in the sense that 

V<£t = 47rVA~ 1 er 4 = -47rVA _1 V • (p + o)u = -47r(p + o)u. 
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After higher order terms being absorbed we obtain 

pi 



2 '' <C(£? + £{) + C(£t +i + £° »). 



Considering 9 t (5.21) and d t (5.22), each term in the RHS cither has been esti- 
mated or can be dealt with in the same manner as the previous cases, and therefore 
the estimates on d 2 d x a and d 2 d x u follows: 

£?■' < C{£f + £f + £l) + C(Ef + s + sl + , + 



Now we move onto two spatial derivative terms, the case i = 2. The only but 
important difference is to use a different weight to close the estimate. It explains 
why the total energy is designed with having different weights according to the 
number of spatial derivatives. Compute d x (5.10) and <9 X (5.11) to get 

(5.23) 

V8 x a = -l d *(P° ±p. { Ut - ^l p -l Vpoa + W<j>+(u- V)« - hV Po } 
w (po + v)~ 75 

_ 7~(Po + v)Hd x u t - ^-p ^Vp d x a + Vd x (p + (d x u ■ V)u + (u ■ V)d x u 

47T (O 

167T -2 

- d x hV Po - -—d x (p 5 Vpo)c - hVd x po} 
75 

V ■ 8 x u = d ; (P ° + r 2 } Wt + Vpo-u + Va-u} 
(5.24) | {P0 + (T) 

{d x a t + Vpo • d x u + Vd x a ■ u + Va ■ d x u + \7d x p ■ u}. 

Po + o- 

~ 2.2 

In order to get right exponents — | + I, | + I of Vd x a and V • d x u in £\ we 
multiply (5.23) and (5.24) by (po + cr) _ i + 5 and (po + cr)ra+5 respectively, and 
square them. Notice that our chosen weight functions are of polynomial type due 
to the behavior of po: (po + a) <~ po = 0(r~ 5 ). Thus as one takes the spatial 
derivative, one gets \d x p \ = 0(r~ 6 ) <~ (p + a)i . So we get the following: 
(5.25) 

J {p Q + a)-i +l \Vd x a\ 2 dx + j (p + a)i +l \ V • d x u\ 2 dx 

<C [ ((po + a)i +l + | dxa ^ \ 2 (po + «) 1+l )(\u t \ 2 + |p -*Vpo| V + |V^| 2 
J (po + a) i" 

+|(u- V)u| 2 + iVpol 2 ^ 2 )^ 
+C [ ((po + a)--* +l + | dx ° A 2 (po + a)-i +l )(\a t \ 2 + |V Po ■ u\ 2 + |V<7 • u\ 2 )dx 

J (po + O) io 

+C J (po + a)i +l (\d x u t \ 2 + |po f Vpo| 2 |^a| 2 + IV^ 2 + \(d x u ■ \7)u\ 2 

+ \(u ■ V)d x u\ 2 + \V Po \ 2 \d x h\ 2 + \d x (po f Vp )| V + \Vd xPo \ 2 h 2 )dx 
+C J (po + <j)-r,+ l (\d x (j t \ 2 + \\7p d x u\ 2 + \\7d x <ru\ 2 + \\7ad x ■ u\ 2 + \\7d x p u\ 2 )dx 
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The first and second integrals are exactly same as the RHS of (5.12). The ones in 
the third and fourth integrals except the potential term have been already estimated 
at the previous steps since each of them contains only one spatial derivative with the 
right exponent of the corresponding weight. The potential part does not produce 
any further difficulty, indeed it behaves better both in weights and in derivatives. 
Notice that Ad x <fr — 4:ird x cr and J R3 d x adx = 0. Thus we can do the same as we 
did in (5.17), (5.18), (5.19) and we get the following estimate similar to (5.20): 

(5.26) J (po + a)i+ l \S/d x 4>\ 2 dx < CS\ + » 

Consequently, we get the desired estimates for d x a and d x u: 

Si 2,2 <C(£f + £l) + C(£l + s + £? +3 _) 

1 5 1 5 

Consider 9 t (5.23) and 9 t (5.24). Taking d t does not destroy the structure of 
equations. As going along the same track, the desired result on d t d x a and d t d x u is 
easily obtained. 

Lastly, for three full spatial derivative terms, namely the case i = 3, j = 3, 
compute <9a;(5.23) and <9 X (5.24). Since we are dealing with one more spatial deriv- 
ative, we have to modify the weights again. Multiply them by (po + er)~To + 2 an( j 
(pa+o-)i + ^ , square and integrate them. Then most terms have been already treated 
before. As for the potential part, noting that Ad x (j) — ^d x a and J R3 d x adx = 0, 
we get the similar estimate to (5.26): 

/ (Pa + a)i +l \Vdl4>\ 2 dx < C£f + * 

At last this finishes the lemma. □ 

Lemma 5.2 shows that any spatial and mixed derivative terms can be estimated 
in terms of time derivative terms with suitable weights, i.e. £i and £i are more or 
less equivalent measurements. Now we take time derivatives which do not destroy 
the structure of the system much and do the energy estimates. 



6. Weighted Nonlinear Energy Estimates 

In this section we perform the nonlinear energy estimates with the utilization of 
a family of symmetrizers of the system. Energy estimates with weights on d\ a and 
df u for < j < 3 are carried out to derive the following key estimate so as to build 
the bootstrap argument which will be discussed in the next section. Throughout 
this section, (5.1) and (5.5) are assumed. 

Proposition 6.1. Suppose £_e + £_z + £_s < 9\ for <t <T where B\ <C 1 is 

5 5 5 

sufficiently small. 

(1) Let I = 0. Then, for any fixed small r\ > 0, there exist C, C v > such that 
~£o < (Cyfc + Vfo + C n (£ )H£^ + C v (£° + £'+£%). 
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(2) Let I < 0. Then, for any fixed small r\ > 0, there exist C, C v > such that 
1 d 



2 dt 



£i < {Cy/Bl + rfiEt + C£ l+ 3 +C v (£ l )i(£ l+{ ^ 2l) )i + C v (£? + £} +£?). 



In particular, if I < — |, then we have 
1 d 



2dt 



Si < (C^Oi + n)Ei + C£ l+ ^ + C^f + C v (£? + e\ + £?). 



Proposition 6.1 will be proven by a series of lemmas in which we will derive the 
estimate on £\ for each j ; each lemma has its own significance and we will need 
all of them to prove the bootstrap argument. Let us start with the simplest case 
j = 0: the zeroth order estimate. 



Lemma 6.2. (£\ P) Let I < 0. There exists a constant C > such that 



Proof. Consider 



Sli 5 jT\-(")dx = 0: 



5.3/ \U 



J (po + ^) ^ + Wcfe + ^ J(po + o-) ■ (po + a)u]ndx 

+ (p + o-) 1+l u t ■ udx + (p + o-) 1+l (u ■ V)u • licte 

(6 ' 1} ir /■ _ 4 o 

+ (po + o-) 1+l [{p + o-) iWa--p 5 Vp o- + h(a,p )]-udx 

+ J {po + o-) l+l \7<j)-udx 
We compute the first three terms in turn; call them (I), (II) and (III). 

<" -55^/0* + " 53(-S + /<» + 

(-") = - tt(~7 + / (Po + (T)-5 + 'V(po + <t) • [(po + 

(6.2) 5 7 4 

- — (po + °~)~ 7 ' +l (po + • Vo-dx 
(HI) =\jj (po + <7) 1+ V| 2 ^ - ^(1 + /(Po + ( r)'a t | U | 2 dx 
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After cancellation (6.1) becomes 
(6.3) 

lj t ^j(P° + °r l+l v 2 dx + j(p +o) 1+l u 2 dx] 

= 1 -^(-Ul)[ (po + *)-i«-2-o*dx + hi + l) ! (po + *f +l -^\u\ 2 dx 
2 15 5 J po + a I J p + a 

- J {po + (r) 1+l (u • V)u • udx - J(Po + ct) 1+ 'V0 • udx 

+ / (Po + a) 1+ '[(- - + l)(po + cr)-^W(po + cr)a + -p 5 Vp <J - h(a, p )] ■ udx 

Next we estimate the potential part. By the Cauchy-Schwartz inequality and 
(5.20), we have 

(6.4) 

J (po + cr) 1+l V0 -udx<^J (po + <7) 1+l -rz\V<j)\ 2 dx + 1 J ' (p Q + a) 1+l +* \u\ 2 dx 
< C£° +JL 

'T 10 

The last term in (6.3) is rewritten as 
(6.5) -|| J(p + a) 1+l [l Po*Vp a + h(a,p )} -udx 

~ — — — 

where h is higher order term including a and Vcr. Recall p 5 p'o — — 5r(l + r ) ~ 

(po + a)~s. When I — 0, we only have cubic terms left in the above and hence we 

are done. The quadratic term when I ^ can be treated as the following: 

(po + cr) 1+ '[p " Vp cr] • udx 

\ j(po + <?) l+l+ ^\u\ 2 dx+ l - J(po + a) 1 + l -^(p- i p' fa 2 dx 



(6.6) 



< 



\ J(po + o) 1+l+ ^u 2 dx + C J(po + a)-i+ l +^a 2 dx 
< C£?, s 



< 
~ 2 



Apply Lemma 5.1 to the first three integrals in the RHS of (6.3). With (6.4) 
and (6.6) the wanted result follows. □ 

For higher order terms, the spirit of details is same as before but we have extra 
terms to deal with. While doing higher derivatives, the necessity of the cooperation 
with mixed and spatial estimates occurs. 

Let us compute 



/*<)•<:) 



dx = : 
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(6.7) 

1 d r 47T 

2~dt 



J (po + a)-i +l (dta) 2 dx + J ( Po + a) 1+l \dlu\ 2 dx} 

I|(-i + 0/(Po + .)-t + ^(^ 

+ ^ l\po + *) 1+l ^\dtufd X 
2 J po + cr 

~%J(Po + vr^diV ■ [(po + a)u}diadx 

- j(po+ v) 1+l d J t [(u • V)u] • d 3 t udx 

-%J(Po + °) 1+l {dtt(Po + a)-iVa] - ^Vpo^a + d^Po)} ■ d\udx 

- j(po+<j) 1+l Vd J t (j)-d J t udx 



For computational convenience, we separate d{ terms from lower derivative terms 
in the RHS of (6.7). Some terms contain unfavorably (j + l) th derivative terms. The 
worst terms seem to come from the third and fifth integrals: — y| J(po + ct)~ s+' V • 

[(p + a)diu]diadx and -*| J(p + a) 1+l [(p + Vd{a - fp 1 V po&{a] ■ d\udx. 
Use the integration by parts to get some nice cancellation: 



(6.8) 

" % I iP0 + CT )" H ' V ' ttP° + ^M^adx 

~ J (po + a) 1+l [(po + ayiVdia - po^a] ■ dfudx 

J (~ + 0(A) + °y^ +l V(po + a) ■ djudladx + (p + a)i +l ^u ■ Vd{adx 



Air 
15 



- (j) 1+l [(po + a)-ivd 3 t a - *V podia] ■ &>udx 



= ^J (po+v) 1+l [I Po § V podia + h(a, Va, p, 8P t a)] ■ & t udx 



NONLINEAR INSTABILITY IN GRAVITATIONAL EULER-POISSON SYSTEM FOR 7 = f 29 

where h{a,Va, p,d{a) = (-| + Z){(Po + cr)~^Vcr + [(p + cr)~? - ']Vp }3 t V. 
Notice that we have the above quadratic term only when I ^ 0. Taking into account 
(6.8) and grouping by similarity, we rewrite (6.7) as following: 



~{g/(A) + ^)-* +, (^) 2 ^ + /(A) + ^) 1+, |^«| 2 ^} 
• 147T, 4 . , x , N _ 4+; a t , 2 



^2 15(-5 +i V^ +ff )"^^^ 



+ ^ /(P0 + a) 1+i -^-|9M 2 dx} 

2 ,/ P0 + 

+{--^ /" (A) + ^)" I+I [^V • u + Vd 3 t a ■ u]d J t adx 

- j (po + [{d 3 t u • V)« + (u ■ V)d 3 t u] ■ &>udx} 

Air f _2 ~ 

(6 9) +{ 15 J (Po + ^l 1 Po^Podfa + h(a,Va,p,d}a) + d}h(a,p )}-dludx} 

+{- J (po + a) 1+l Vd{<P ■ Oludx} 

H-% j (Po + ]TV. [di-\p + a)diu]diadx 

-I(po+ <y) l+l £p*u • v)3T«] • diudx 

i=l 

-%J(Po + £0T(Po + «)-^dlV*} ■ diudx} 
=(/) + (77) + (7/7) + (IV) + (V) 



Note that (V) is alive only for j = 2, or 3. First three groups have the exactly 
same structure as (6.3). Each term can be easily estimated as in Lemma 6.2. For 
(I), by using Lemma 5.2, we get immediately 



(6.10) 



(I) < Cy/elel for all j. 
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(6.11) 



Other (j + l) th derivative terms are in (//). After integrating by parts, the third 
and fourth integrals become 

~ if j {fi0 + fT ) _l+ '[^ f7V ' u + Vd l a ■ u]d J t adx 

27T. 4 f V(po + O") U , 2 , 

+ tt(-- + / (po + 5+ ' K .J_ ■ -rWofda, 

ID D J (p + a) io (p + er) io 

J (p + a) 1+l [(diu • V)u + (« • V)$u] • S^'uda; 
= - J( Po + cr) 1+l (dlu • V)« • ^'uda; + \j(Po + cr) 1+l (V • u)|#'u| 2 da; 
+ -5- / {po + <?) 1+l - H \ .1 ■ - — — -\d J t u\ 2 dx 

1 J (Po + cr) 10 (po + a) 10 

Therefore (6.11) with Lemma 5.1 gives rise to: 
(6.12) (II) < Cy/9l£i for all j. 



(Ill) is similar to (6.5) in the zeroth estimate and so it can be treated in the 
same way. If we do the same as in (6.6) and use Lemma 5.1, we get the following: 
for each j, 

(III) <CJ¥ x Ei when I = 0, and 

(6.13) . ^_ . 

(Ill) <CSf, 3 + Cyf¥ x e 3 , when I ^ 0. 
'+ 10 

The potential part (IV) and cubic terms (V) are somewhat new, complex and 
they'd rather be done separately according to different j's. Before we split the 
cases, from the dynamics of S7d{(j> : Vd{(j> = —4ird{~ 1 [(po + cr)u], we reduce (/V) 
to the following: 

(6.14) (IV) = 4ttJ (po + a) 1+L d{u ■ d^Kpo + <?)u]dx 

Now let j = 1. Here is the first order d t estimate. 
Lemma 6.3. (£}) For any small r\ > 0, there exist constants C, C v > such that 

\j t £l < (C^Ol + v)£o + C& f° r 1 = 0, 

< (Cjel + n)£l + C£l + s_ + C v £? +2 for I < 0. 

Proof. Since (V) has no effect on j = 1, it is sufficient to take care of (IV). The 
potential part is shown to be even better in terms of derivatives as we can expect 
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in (6.14). By the Cauchy-Schwartz inequality, we get 
(IV) = 4tt J {po + <jf +l d t u ■ udx 

(6 ' 15) <V f(po + a) 1+l \dtu\ 2 dx + C n J( P0 + af+ l \u\ 2 dx 

< r\£\ + C n £i +2 , for any small 77 > 0. 
Thus (6.10), (6.12), (6.13) and (6.15) give the desired result. □ 

The second order d 2 estimate can be done in the same spirit. We have extra 
cubic terms to deal with from (V). Lemma 5.2 plays an important role. Let j = 2. 



Lemma 6.4. (£f) For any small w > 0, there exist constants C,C V > such that 

1 d 1 

2 Jt £ I <(cVo~i+ v)£% + Cr,J2 4 V 1 = 0, 



i=0 

1 



\j£f < {Cy/fx + n)£ 2 + C£ 2 + s_ o +a n J2 Z\ for I < 0. 

i=0 

Proof. The potential part (6.14) can be computed like (6.15). By Lemma 5.1, 

(IV) = 4tt ( (po + o-) 2+l d 2 u ■ \d t u + -^— u ]dx 
(6.16) J po + o- 

< r\£ 2 + C v (£l +2 + ^1^+2)1 fo r an Y small 77 > 0. 

Terms in (V) for j ~ 2 are at least cubic including mixed derivatives. We can take 
the sup for the lowest, first derivative term and then we end up with manageable 
quadratic terms. In order to see how it works, we illustrate the estimate on the 
first term in (V): 

(6.17) 

J {Po + o-)-i +l V • [dtad t u]d 2 adx 

= f(p + o-)r,+ l dta (V • d t u)d 2 adx + [ (p + <j)~^+l 9tU . Vd t ad 2 crdx 
J (Po + o-) J (po + o-)"' 

<C^9[{ J (po + <j)t +/ |V • d t u\ 2 dx + J(po + a)-i+ l (d 2 a) 2 dx} 

+ C^F 1 { J (po + a)-l+ l \Vd t a\ 2 dx + J (p + a)-i +l (d 2 a) 2 dx} 
<Cy/6l{£i* + £?) 

We have used the Cauchy-Schwartz inequality at the first inequality. Notice the 
changes of the exponents in weights. We can do the same to other terms in (V) . 
Ultimately, applying Lemma 5.2, we get the following: 



(6.18) (V)<C^0l£ 2 +cY,£ l+ 



1 

k 

5 



fc=0 
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Note that £i +k < C£i for k > 0. Thus (6.10), (6.12), (6.13), (6.16) and (6.18) 
give the desired result. □ 

To finish the proof of Proposition 6.1, only j — 3 i.e. df case is left. The difficulty 
is to handle the weighted new cubic terms in (V") of which each factor is at least 
second derivative of a and u and hence we cannot utilize Lemma 5.1 directly. To 
overcome it, we introduce the weighted Gagliard-Nirenberg inequality. This job is 
done in the next lemma. One can see that (IV) and other terms in (V) for j ' = 3 
can be treated similarly as in (6.16) and (6.17). Therefore, the following lemma 
finally establishes Proposition 6.1. 

Lemma 6.5. Let I < 0. For any small fixed -q > 0, there exists C v > such that 

J (p + a)-^ +l \d^(rdtV ■ udfa\dx, J (p + a)-^ +l \d t Va ■ d?ud?a\dx 
J (Po + <j)-r,+ l \d 2 t ad t V<T ■ dfu\dx, J (p + cr) 1+l \d 2 u ■ Vd t ud?u\dx 

are bounded by n£f + C v (£f)i (£f + (_ 3 _ 2 ;))^ + ^vi^f+s. ) 5 (^i+(-3-2i))' ■ ^ n partic- 
ular, if I < — |, they are bounded by r\£i + C v (£i) 2 . 

Proof. First we split each term in df term and d 2 term by the Cauchy-Schwartz 
inequality. In order to take care of L 4 norm of d 2 terms we use the Gagliard- 

3 1 

Nirenberg inequality ||/||l4 (r3) < allWII^^)!!/!!^^)- Since the inequality com- 
plies well with the localization and the weight functions are nice, using a partition 
of unity, in our case we get the weighed version of the Gagliard-Nirenberg inequality 



/>oc />oo />oo 

/ w k f 4 r 2 dr<C( w^\Wf\ 2 r 2 dr)i( w%f 2 r 2 dr)^ 

(6-19) Jo h J » 

+C( w^'f 2 r 2 dr)^( w(f 2 r 2 dr)^ 
Jo Jo 

where Wk = (1 + r 2 ) 1 ^ ~ (p + cr) 1- ^ and |a + \fl = |a' + \f3' = 1. Its proof is 
given at the end of the argument. Only the first and the last term are treated in 
this proof. The other cases can be estimated in the same way. 



(6.20) 

(p + a)-r,+ l \d 2 ad t V ■ udfa\r 2 dr 



<r) (p + o-)-i+ l (dfo-) 2 r 2 dr + C n (p + a)-^ +l \d 2 ad t V ■ u\ 2 r 2 dr 
Jo Jo 

r-OO 

/ (p + o-)-t +l \d?ad t V -u\ 2 r 2 dr 
Jo 

<o / (po + <y)-^ +l (d 2 nfr 2 dr+- / (p + a)i +l (d t V ■ u) 4 r 2 dr 
2 Jo * Jo 
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(6.21) 

pOO 

/ (po + <7)-T+'(afc)Vdr 

JO 



poo poo 

<C{ / (p + a)-i +l \d 2 Va\ 2 r 2 dr}i{ / (p + a)-™~ 1 (d 2 a) 2 r 2 dr}i 
Jo Jo 

poo poo 

+ C{ / (po + a)-*+'(^)Vdr}i{ / (po + a)-^- l (d 2 a) 2 r 2 dr} } . 



JO 



<C(£t)H£f +( _ 3 _ 2l) )->+C(£f +1 JH£f +{ _ 3 _ 2l) ^ 
(po + ff) l+ '|a t V-w|Vdr 

<C{ (p + a)i+ l \d t V(V -u)\ 2 r 2 dr}i{ (p + a)-i~ l \d t W • u\ 2 r 2 dr}* 
Jo Jo 

poo poo 

+ C{ (p + a)i +l \d t V-u\ 2 r 2 dr}i{ (p Q + a)-i~ l \d t V ■ u\ 2 r 2 dr}i 
Jo Jo 

Here is the verification of each exponent: — + Z = |(— | + Z) + ^(— ^ — Z) and 
-f -Z = -|+Z+(-3-2Z); |+Z = §(g+Z) + ± (-|-Z) and -|-Z - |+Z+(-3-2Z). 
Note that — — Z > — | + Z for Z < — §. Now let us look at the last term. We go 
through the similar computation as in (6.20) and (6.21). 

(6.22) 

pOO p 00 

/ (p + a) 1+l \d 2 u-Vd t udfu\r 2 dr < V / (p + a) 1+l \dfu\ 2 r 2 dr 
Jo Jo 

poo poo 

+ C V { (po + ^ +l \d 2 u\ i r 2 dr+ (p +a)i +l \Vd t u\ 4 r 2 dr} 
Jo Jo 



poo 

/ (p +a)i +l \d 2 u\ 4 r 2 dr 
Jo 



(6.23) 



pOO pOO 

<C{ (p + a)i+ l \d 2 Vu\ 2 r 2 dr}i{ (p + a)- 2 - l \d 2 u\ 2 r 2 dr}i 
Jo Jo 

poo poo 

+ C'{ {p + (7)i +l \d 2 u\ 2 r 2 dr}l{ (p Q + a)- 2 ' 1 \d 2 u\ 2 r 2 dr}^ 
Jo Jo 



<C(ff)t(ff +( _ 3 _ 2I) )* + C(£? H )*(5/V ( _3_ 2I) )* 
{p Q +<j)* +l \Vd t u\ 4 r 2 dr 
<C(£f)i(£f +( _ 3 _ 2;) )nC(£f + ,)i(£ i 2 +( _ 3 _ 20 )^ from (6.21) 



I 



Observe that § + Z = §(§ + Z) + I(-2 - I), § + Z = |(| + Z) + §(-§ - Z) and 
-2- Z = 1 + Z + (-3-2Z), -f -Z = | +Z+ (-3- 21). 

The only missing part is the proof of the weighted Gagliard-Nircnberg inequality. 
Here it comes: 
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Proof of (6.19): We choose a partition of unity {(p n }n>o as following: 

fo(r) 

' 



1 - r if < r < 1 
if r > 1 



<Pn(r) = < 



if < r < 



re— 1 -f re — 1 
2 11 2 



if f < r < 2±i 



«±2_ r if«±i< r <«±2 
if r > «±2 



It is easy to check < <p n < 1, supp = [ IL 2^, n Y^]{ supp </? = [0,1]), 

Erelo ^"( r ) = 1 for a11 T ^ °) and l^nMI < 1 a - e - 

^OC (-00 °° 

/ {l + r 2 ) k fr 2 dr = [ V ^„(r)] 4 (l + r 2 ) fe /Vdr 



(6.24) 



<CY, ^ n (r)(l+r 2 ) k fr 2 dr 
re^ 



We only consider k > 0. Other cases can be proven in the same manner. Note 
that (1 + r 2 ) k < (1 + { 1J 4^) 2 ) k on supp <p n . First we localize the half real line 
according to the partition of unity. Since weights are monotonic, they can be 
localized as well. And then we apply the Gagliardo-Nirenberg inequality. 

(9.1)<C^(l + (^) 2 )^ °^/Vdr 

<^E( 1 + (^) 2 )' £ ( J / o [My] 2 r 2 dr)Hj o <p n 2 f 2 r 2 dr)i 
<^E( 1 + (^) 2 )' £ ( J / [M'f + yjf]r 2 dr)Hj o ^fr*dr)h 

n + 2 n + 2 

n n+2 n + 2 

+ cE( 1 + (^) 2 ) fe ({f + /J }fr 2 dr)H[_ 2 i f 2 r 2 dr)^ 



+i 

2 

+2 



<C]T(y i (l + r 2 ) a f' 2 r 2 dr)i(J ^ (1 + r 2 f f 2 r 2 dr)^ 

n n+2 n+2 

+ CJ2({[\ + [ \ }{l + r 2 ) a ' fr 2 dr)i{( * (1 + r 2 f'f 2 r 2 dr)i 

<C{ (l + r 2 ) a f' 2 r 2 dr)i( (1 + r 2 ffr 2 dr)^ 
Jo Jo 

f-OO f-OQ 

+ C( (l + r 2 ) a ' f 2 r 2 dr)i( (1 + r 2 f f 2 r 2 dr)^ 
Jo Jo 

In the above C is a generic constant. Note that |a + i/3 = §a' + i/?' = fc.D 
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7. Nonlinear Instability 

Now we are ready to prove the bootstrap argument. The proof completely de- 
pends on the estimates in Section 6 and the Gronwall inequality. 

Proposition 7.1. Let v{t) = (^j) be a solution of the Euler-Poisson system (1.9) 
and (1.10). Let 1+ < —3 be given. Assume that 

VSTAQ) < Co<5 and ^/ig(t) < CoSe^ 5 * forO<t<T. 
Then there exist C$, 8 n > such that 

if 0<t< min{T,T 5 }, then \JWm < C b 8e^ < C 5 6> , 
where T s — -j= In % . 

Proof. In Proposition 6.1, Lemma 6.2, 6.3 and 6.4, choose 6\ and rj small enough 
so that C\[Wi + rj < -^y^. Therefore there exist constants C\ > 0, C 2 , C 3 > such 
that for I < 0, 





1 d o 

2 (ft ' 


< 


2 


+ C<2,£,. 3 
10 




ft) 


1 d ! 

2 eft ' 


< ■ 


2 


+ C\E], 3 
i+ 10 


+ C2^i° + 2J 


(q) 


1 d £ 2 

2<ft 1 




2 


+ C\Ef, 3 
10 


+ C 2 (£°+£ ; 1 ), 




Id 
2dt 1 




2 ^ 


f Ci£,,3_ - 

1+ 10 


f C 2 (£ ; ° + £} + £?) + C 3 (£ l )i(£ l+( ^ 2l) )i 



Note that G\ = when Z = 0. Define T* by 



T* = sup{t : £,(«) < mm{6 1: [0,t], Z* < Z < 0}. 

4C/3 

Here 6*i is a small constant coming from the smallness assumption to guarantee 
that the nonlinear energy estimates work. 

Let < t < min{T,T*}. Then since y/£$(t) < CoSe^ 1 by the hypothesis, 
from (a;), we get 



(aj) => \J£i(t) < C'oSe^ 1 for all I < 

by the standard Gronwall inequality. We use Cq as a generic constant. Consider 
the following diagram: 

(bo) \jfi<^£l + C 2 £°^^(t)<C' 5e^ 

Likewise, for all k > 0, we have 

(6_^ fe ) => ^/^(t) < CSfc^«. 
By the same bootstrap argument with (c_3_ k ) where k > 0, we have 
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Now we move onto (di). Let us start with I = 0. Recall that C\ = 0. Since 



(£ ) I (S _ 3 ) 2 < firstly we get 

-£o < V^fo + 2C 2 (£ ° + Si + S 2 ) + 2C 3 (S )HS^ 



4C 3 

d_ 

< M^£ + 2C 2 C 2 S 2 e 2 ^ t 

By the Gronwall inequality, we have the exponential growth with the exponent 
on Sq. As the previous cases, form a sequence {(d_3_ k )}k>o and use the 
Gronwall inequality to get the following 

(7.1) Si(t) < Cl8 2 e 2 ^\ for < t < min{T,T*} 

for any I < and some constant C4. And in success by Lemma 5.2, for any I < 0, 
we also have 

(7.2) Si(t) < Cl8 2 e 2 ^\ for < t < min{T,T*}. 



Now choose 0q such that (C^Oq) 2 < min{0i, j^}- We consider following two 
cases: 

(i) T 5 < min{T, T*}; in this case, the conclusion follows without any extra work. 

(ii) T s > min{T,T*}; then T < T* < T 5 . If this is true, then again the 
conclusion is trivial. We show this is the only possibility. If not, we have 

T* <T < T s . Letting t = T*, from (7.2) and the definition of T s , we get 

Si(T*) < C 2 5 5 2 e 2 ^ T * < (C 5 Sexp^ TS ) 2 = (C 5 9 ) 2 . 

But this is impossible by the choice of 9o since it would contradict the definition of 
T*. This establishes the proposition. □ 

From now on we regard S > as an arbitrary small parameter and 9 as a 
small but fixed positive constant independent of S. Recall that T s is defined by 
9 = (5expv / ^ T<5 or equivalently T s = -^== In f . 

Proof of Theorem 1.1. Let vq = (^°) be a growing mode for the linearized 
Euler-Poisson system (2.1) and (2.2) obtained in Section 2. Normalize vq such that 

INlN^T 2 (Mvo + ll^ ) = i- 

By the behavior of ^0 in Proposition 3.4 and the relation (2.4) between 0o and "00; 
we know \Svo\ is relatively small compared to p a ; p + S(f>o ~ 7Po where 7 is close 
to 1, and / °° <f> r 2 dr = 0. We may assume ||po < Po + <50o < %Po- Let 



/■oo ^ />oo 

4tt / _( (0o + ( 5 ( /, )-|+ / ^r 2 rfr + 47r / (p 
Jo to Jo 



+ J^o) 1+ Vo^dr = a 2 < 00. 

Now solve the Euler-Poisson system with a family of initial data v\t=o = Svq. 
The continuity equation gives rise to <r s r 2 dr — 0. Denote the corresponding 
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£;-solution by v(t) = v s (t) = (^|*|)- It can be written as 

v(t) = Se^uo + f C{t- t)N{j)6t. 
Jo 

where C is the solution operator for the linearized Euler-Poisson system and N is 
nonlinear part. 

— £ 



N = 



where h = (p + a) 5 - p 5 + |p 5 a. 
Define T by 

T = sup{s : for < t < s, < max{3, ajSe^ 1 "*}. 

Then by Proposition 7.1, there exist Ci and f? n > such that 
for < t < min{T,T 5 }, 

y/ei(t) < \fii(t) < dSeyP 5 * < Ctfo. 

Note for sufficiently small 9 0l by Lemma 5.1, it means that | p ° +IJ i \ « 1 i.e. 
po + <J S behaves like po- So we can find a small constant (i, < (i < | such that 

(? 3) (1 - /?) 2 Po~ f < (Po + ^)"* < (1 + /?) 2 p ~ f 

(1 - /3) 2 po < Po + ^ < (l + /3) 2 Po- 

To emphasize which functions wc deal with, we denote £® being plugged / = (jj 
in but its weight part unchanged by ||/||y: 

ll/ll * r = 47r / 7^o + ^r l (.A)V rfr + 47r j / o (P0 + ^)(/ 2 )Vdr 
In this notation, by (7.3), for < t < min{T,T s }, 
(7.4) ||(5 e v^Vo||Y > (1 - /3)<5e^'||„ || = (1 - /3)£e^*. 

For the nonlinear parts, from the linearized estimates in the Lemma 4.1, 4.2 and 
4.3, for < t < min{T, T s }, we have 

(7.5) 

||^(t)-fcv^V ||Y=|| / C(t-T)N(T)dr\\ Y 
Jo 

<C f e^ (t - T H\\N(r)\\ Y + \\d t N(r)\\ Y + \\d 2 t N {r)\\ Y )d T 
Jo 

<C /' e VP.(«-r) { |_g_| + |_gj_| + | Va 

Jo PO + cr po + cr (p + (j)To 



+ L " ,J + I + |VU| } WSl+S^dT 

(Po + cr) 10 (p + cr) 10 

t 



<C f e^ {t - T \5e^ T ){5e^ T )dT 
Jo 

<C 2 (5eV^*) 2 
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where C2 is a constant. At the second inequality wc have used Lemma 5.2. The 
next inequality follows from Proposition 7.1 and Lemma 5.1. 
Now if necessary, fix 9 sufficiently small such that C 2 < 

Claim. T s < T. 

Proof. If not i.e. T s > T, by (7.4) and (7.5) 

\W 5 \\y(T) < ||<5e^V ||y(T) + \\u 5 - Se^isoMT) 

< (1 + 0)5e^ T \\is o \\ o + C 2 9 Se^ T 

< ?±l S e^ T < 2Se^ T 
which would contradict the definition of T. 

Once we have T s < T, again by (7.4) and (7.5), 

V^(T 5 ) > (1 - (3)5e^ T& - ^lse^ T * = ^-j^0 > 0. 
Set 6 = ^-^-6q. This finishes the proof of the theorem. □ 
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